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We explore the precision with which the Einstein Telescope (ET) will be able to measure the 
parameters of intermediate-mass-ratio inspirals (IMRIs), i.e., the inspirals of stellar mass compact 
objects into intermediate-mass black holes (IMBHs). We calculate the parameter estimation errors 
using the Fisher Matrix formalism and present results of Monte Carlo simulations of these errors 
over choices for the extrinsic parameters of the source. These results are obtained using two different 
models for the gravitational waveform which were introduced in paper I of this series. These two 
waveform models include the inspiral, merger and ringdown phases in a consistent way. One of the 
models, based on the transition scheme of Ori & Thorne 0], is valid for IMBHs of arbitrary spin, 
whereas the second model, based on the Effective One Body (EOB) approach, has been developed 
to cross-check our results in the non-spinning limit. In paper I of this series, we demonstrated 
CNj ■ the excellent agreement in both phase and amplitude between these two models for non-spinning 

black holes, and that their predictions for signal-to-noise ratios (SNRs) are consistent to within ten 
^ | percent. We now use these waveform models to estimate parameter estimation errors for binary sys- 

|"t | ! terns with masses 1.4M +1OOM , 1OM +1OOM , 1.4M +5OOM and 1OM +5OOM and various 

choices for the spin of the central IMBH. Assuming a detector network of three ETs, the analysis 
shows that for a 1OM compact object (CO) inspiralling into a 1OOM IMBH with spin q — 0.3, 
detected with an SNR of 30, we should be able to determine the CO and IMBH masses, and the 
IMBH spin magnitude to fractional accuracies of ~ 10~ 3 , ~ 10~ 3 ' 5 and ~ 10~ 3 , respectively. We 
' also expect to determine the location of the source in the sky and the luminosity distance to within 

Qi ~ 0.003 steradians and ~ 10%, respectively. We also compute results for several different possi- 

ble configurations of the detector network to assess how the precision of parameter determination 
^q' depends on the network configuration. 
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"sj" ■ I. INTRODUCTION 

O , 

Significant progress has been made over the last few years in the quest for the first direct detection of gravitational 
waves (GWs). During the most recent published data-taking run of the ground-based GW detectors five instruments 
were operational. The three LIGO detectors started their Science Run 5 (S5) in November 2005 0|. Subsequently, 
the GEO 600 detector Q joined the S5 run in January 2006 and the Virgo detector [3] started its Science Run 1 
(VSR1) in May 2007. All five detectors worked in triple-coincidence until the beginning of October 2007. The data 
collected during the rest of the S5/VSR1 run, referred to as "S5y2/VSR1", was the first long-term observation with 
the worldwide network of interferometric detectors. Towards the end of 2010, the LIGO and Virgo detectors were 
taken offline and the upgrade to the Advanced detectors began. If the design sensitivities are achieved, Advanced 
LIGO and Virgo will increase the detection range of future searches by as much as a factor of 10 over those reached 
in S5y2/VSR1, so that the monitored volume of the universe will increase by a factor of ~ 1000. At this level of 
sensitivity, the detection of GW signals from binary mergers should become routine. 

Seismic noise will still limit the sensitivity of the advanced detector network LIGO-Virgo-GEO to frequencies J> 10 
Hz || . Hence, upgraded ground-based advanced detectors will not be very sensitive to signals from coalescing 
binaries in the ~ 10M Q — 1000M© range, which generate GWs predominantly in the 0.1-10 Hz frequency band. 
These events will also be inaccessible to the space-based GW detector LISA, as this instrument will be most sensitive 
to events with total mass ;> 1000M Q , i.e., signals in the O.lmHz — 0.1Hz band. Studying GW sources in 0.1-10 Hz 
frequency band is needed to address an important, yet elusive, problem in astronomy, that of the existence of black 
holes of intermediate mass (IMBHs). There has been a long standing suspicion that these objects may form in the 
centres of dense stellar clusters 0-111 • Some evidence supporting this possibility has come to light over the last decade, 
mainly from X-ray and optical observations 10]. Detection of these IMBHs, particularly in dense stellar clusters, 
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may be important as they would have a direct impact in the dynamical evolution of the host cluster and lead to the 
generation of GWs. 

There are two types of data that suggest the existence of IMBHs: i) ultraluminous X-ray sources have been observed 
in galaxies that are not associated with AGN, but which have fluxes which would imply accretion rates many times 
the Eddington limit, if the accreting object had a mass M <; 20M Q ; ii) the stellar kinematics in the centres of 
several globular clusters, e.g., M15 in the Milky Way, and Gl in M31, show evidence for an excess of dark mass in 
their centres. There are two distinct channels through which IMBHs might have formed. Firstly, they could form in 
the early universe through the collapse of massive, low-metallicity "Population III" stars. In the cold dark matter 
framework, galaxies formed through a series of mergers between initially low-mass objects that condensed at early 
cosmic times. The black holes seen today similarly grew through a combination of mergers and accretion, from initial 
seeds. In the Pop-Ill star model, the initial seed black holes are of intermediate mass, IOOMq — 1000M Q , but the 
initial seeds could also have been 'heavy', ~ 10 5 M Q . The discovery of IMBHs present at early cosmic times would 
be a strong discriminator between these alternative models. The second proposed channel for IMBH formation is the 
runaway collision of stars in the centre of a dense stellar cluster [ill Il2j | . If IMBHs are found to be common in globular 
clusters, this would be an important probe of the processes taking place in such dense stellar systems. Observations 
in the electromagnetic spectrum will improve in coming years, but it is likely that the first robust mass determination, 
and hence the first convincing proof of the existence of IMBHs, will come from GW observations [lCJ. 

In order to detect GWs from systems containing IMBHs, we require a detector operating in the 1 — 10Hz frequency 
range. A design study is currently ongoing within Europe for a third-generation ground-based GW detector, the 
Einstein Telescope (ET) pi. [T3|. to follow on from the Advanced detectors. The target is a sensitivity ten times 
better than that of the Advanced interferometers, and a frequency sensitivity band that stretches down into the 
1 — 10Hz range. The ET will be an outstanding tool to address problems in fundamental physics, cosmology and 
astrophysics. It will have the capability to do the same type of science as the Advanced detectors, but with better 
sensitivity, and to do new science through the detection of sources at lower frequencies. The IMBH sources fall 
into this second category, and in this paper we will consider GWs generated from intermediate-mass-ratio inspirals 
(IMPJs), which are the mergers of IMBHs with lower mass compact objects, namely neutron stars or stellar-mass 
black holes. 

For ET to detect GWs generated during IMRIs, we need accurate waveform models in the intermediate-mass-ratio 
regime, which are not yet available. Source modelling is an ongoing effort that has reached a fairly advanced stage 
for comparable-mass-ratio and extreme-mass-ratio systems. In the comparable-mass regime post-Newtonian (PN) 
theory has provided waveform templates for the inspiral phase evolution, while numerical relativity (NR) can now be 
used to model the final few cycles of inspiral and the plunge and merger [l5j | . At the other extreme we have extreme- 
mass-ratio inspirals (EMRIs) with mass ratios ~ 10 _6 -10 -4 . In this regime, the smaller object may complete several 
hundred thousand orbits at a velocity that is a significant fraction of the speed of light before crossing the horizon of 
the central massive black hole. Modelling the gravitational radiation of EMRIs using NR is presently impractical due 
to the excessive computational cost of simulating so many orbits. However, we can model these systems accurately 
using black hole perturbation theory, treating the mass-ratio as a small expansion parameter [l6l - ll9j | . The sources we 
shall now examine lie somewhere between these two regimes. In paper I of this series [20| we developed two waveform 
models that include the inspiral, merger and ringdown phases using the best of what is currently available. In the 
first model, valid for CO inspirals into IMBHs of arbitrary spin, we used the transition scheme developed by Ori 
& Thorne ll|. In order to build up confidence in this waveform model, we also used the effective-one-body (EOB) 
approach [U [52] to construct a different waveform model for non-spinning systems. The two models showed an 
excellent level of agreement in both the phase and the amplitude. We used these models to compute signal-to-noise 
ratios (SNRs) for a selection of IMRI systems. We compared results from the two models for non-spinning systems 
and found that the predictions were consistent to better than ten percent. We used these SNRs to estimate the number 
of events that could be detected by ET and found that this could be as many as hundreds of IMRI coalescences per 
year up to redshift z ~ 5, if the systems are predominantly light and with high spin, (i.e., an IMRI with redshifted 
masses, [10 + 100]M Q , and with IMBH spin parameter q = 0.9), or as many as tens of IMRI coalescences at redshift 
z ~ 0.3, if the systems are predominantly heavy and with low spin ([1.4 + 500]A/ Q , with q = 0) [2Cj. In this paper, 
we shall address a complementary problem, namely, to calculate the precision with which an ET network will be able 
to estimate parameters of any IMRI sources that are detected. This is an important problem because an accurate 
determination of the masses of a merging binary and a measurement of the luminosity distance at which the merger 
is taking place will be useful for extracting science from the observations. Since the binary systems are typically 
short-lived, a single ET will not be able to localise a source on the sky, so we will assume the existence of a detector 
network. We will consider several different network configurations in order to assess how this affects the science, which 
will be useful input for future design decisions. 

This paper is organized as follows. In Section [II] we outline the assumptions we shall adopt for the detector and 
detector networks, and describe the binary systems to be used in the subsequent analysis. In Section |TTT] we present a 
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brief description of the waveform models, introduced in paper I, that will be used again here to estimate parameter 
estimation accuracies. In Section HVl we summarise some basics of signal analysis that will be relevant for our studies 
and briefly describe the noise model we shall use for our studies, "ET B" [23|. In Section [V] we will use our waveform 
models, as described in paper I 20], to estimate the parameter errors that would arise due to detector noise, using 
the Fisher Matrix formalism. We present results for a 3-ET network for 12 sample IMRI systems and use four of 
these to explore how the results change under alternative configurations of the network. Finally, we summarize our 
work and discuss the implications of these results in Section [VII 



II. ASSUMPTIONS 



A. Einstein Telescope Design 



We follow the same assumptions about the ET sensitivity that we made in paper I. The design sensitivity is usually 
quoted for a single right-angle interferometer with 10km arms. In this paper we shall use the "ET B" sensitivity 
curve, as described in [23j . which was the most up-to-date sensitivity curve available when this work started. The 
corresponding amplitude spectrum is plotted in Figure [TJ 
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FIG. 1: Sensitivity curve for the Einstein Telescope. The Advanced LIGO noise curve is also shown for reference. 



It is expected that ET will operate with an improved sensitivity over advanced detectors in the frequency range 
l-10Hz. Assuming a low frequency cut-off of 1Hz might be too optimistic as it is a significant technological challenge 
to push the sensitivity toward lower frequencies. ET may therefore only achieve sensitivity down to ~ 3Hz or higher. 
Given this uncertainty, we adopt a conservative approach and use a low frequency cut-off of 5Hz. 

We take the response of a "single ET" to be that of two right-angle interferometers, coplanar and colocated but 
offset from one another by 45°. The currently favoured ET design is for a triangular configuration comprising three 
10km detectors with 60° opening angles in a single facility, as this has lower infrastructure costs [H|,[l4j]. The response 
of the triangular configuration contains the same information as that of the two right-angle interferometers (assuming 
uncorrelated noise), but its sensitivity is a factor of 3/(2-\/2) ~ 1.06 higher. As in paper I, we ignore this factor since 
it is small compared to other uncertainties in the design. 

IMRI signals are short-lived, so to enable extrinsic parameter estimation we need to assume the existence of a 
network of detectors. Our optimistic reference case will be a network of three detectors, each with the sensitivity of 
a single ET. We assume the detectors are sited at the geographic locations of Virgo, LIGO Livingston and in Perth 
(Australia). This configuration is highly optimistic, so we shall also explore how parameter estimation accuracies are 
modified for four additional network configurations, C1-C4. These configurations are, CI: one ET at the geographic 
location of Virgo; C2: as configuration CI plus a right-angle detector at the location of LIGO Livingston; C3: as 
configuration CI plus another ET at the location of LIGO Livingston; and C4: as configuration C2 plus another 
right-angle detector in Perth. We will denote the reference 3-ET network as configuration C5. 
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B. Sample IMRI systems 

We will present parameter estimation errors for the same twelve binary systems that we used in paper I [2(3] ■ These 
correspond to all possible combinations between four sets of component masses — 1.4M + 1OOM , 1.4M + 5OOM , 
1OM + 100M Q , and 10M Q + 5OOM — and three different values for the spin of the central IMBH — q = 0, 0.3, 0.9. 
The mass and spin distributions of IMBHs are extremely uncertain, so we choose these systems to cover the range 
of possibilities. The small object masses are chosen to represent mergers with black holes (1OM ) or neutron stars 
(1.4M ). 

As discussed in paper I [20], if IMBHs gain mass through a series of minor mergers with smaller objects, then their 
angular momenta will undergo a damped random walk. Under such an evolution, it is expected that IMBHs in the 
mass range 1OOM — lOOOAf would end up with moderate spin parameters, q ~ 0.3 [2J]. We will regard this as 
the 'typical' case and will use only that spin value when we explore how the parameter estimation accuracies depend 
on the network configuration. For the other spin parameters, we shall only quote results for the optimistic 3 ET 
network and these should be regarded as 'best-case' values for the accuracy with which we may be able to determine 
the parameters. 



III. IMRI WAVEFORM MODELLING 



We will use the two waveform models developed in paper I of this series to explore the precision with which the 
various ET network configurations will be able to determine the source parameters. Full details of the waveforms are 
given in [20| . but we briefly summarise the models here. 

Both waveform models include inspiral, merger and ringdown in a consistent way. We model the inspiral phase 
evolution, in both models, using the "numerical kludge" waveform model described in [l^]. The "numerical kludge" 
approach consists of computing the orbital trajectory that the inspiralling body follows in the Boyer-Lindquist coor- 
dinates of the Kerr spacetime of the central black hole. We assume an adiabatic evolution in which the inspiralling 
body is moving on a geodesic and use the fluxes of energy E and angular momentum L z derived in (25| to evolve the 
parameters of the geodesic. In this paper, we restrict our attention to circular and equatorial inspirals. Subsequently, 
we numerically integrate the geodesic equations along this inspiral trajectory and obtain the Boyer-Lindquist coor- 
dinates of the small object as a function of time. Finally, we construct a gravitational waveform from the inspiral 
trajectory using a flat-spacetime wave-emission formula. This scheme, which was developed for the modelling of 
EMRI systems for LISA, has several useful features: a) the waveforms have been checked against more accurate, 
Teukolsky-based, waveforms for test-particles on geodesic orbits and the overlap exceeds 0.95 over a large portion of 
the parameter spa ce fill ; b) they are computationally inexpensive; c) conservative self-force corrections to this model 
have been derived for Kerr circular-equatorial orbits at 2PN order. Nonetheless, this model is, at present, incom- 
plete in several ways, e.g., conservative corrections are not yet known for generic orbits; the phase space trajectories 
are approximate, although they have been matched to Teukolsky based evolutions; and the waveform is constructed 
from the trajectory using a flat-spacetime wave-emission formula. Nevertheless, this approach should capture the 
main features of the inspiral waveform accurately. 

The two models differ in their treatment of the plunge and merger phases. The first model, valid for IMBHs of 
arbitrary spin, employs the transition scheme developed by Ori & Thorne [l[ to smoothly match the inspiral phase onto 
the transition to plunge. In the spirit of the "numerical kludge waveform" , we derive an approximate waveform from 
this transition trajectory using a flat-spacetime emission formulae applied to the transition orbit and plunge geodesic 
given by the transition model [201 ] . The inspiral waveform includes the dominant modes £ — m = 2, I = — m = 2 only 
and so, for consistency, we included only the same two dominant modes in the expression for the waveform generated 
during the merger phase. 

The second model, used for non-spinning central black holes, uses the EOB prescription [2l|, [22j to describe the 
plunge and merger phases. The EOB model is a framework in which the motion of a binary is represented as the 
effective motion of a test object in a background. Although we use the EOB model only in the non-spinning limit, 
an extension of the EOB scheme does exist which includes leading-order spin-orbit and spin-spin dynamical effects 
of a binary system for an "effective test particle" moving in a Kerr-type metric [26|], and next-to-leading-order 
spin-orbit couplings [13] ■ However, it has been recently found 28] that it is not straightforward to include higher- 
order non-spinning PN couplings, such as the 4PN and 5PN adjustable parameters that were recently calibrated to 
numerical relativity simulations for non-spinning systems [1^, [30] , using these Hamiltonians [26|, HtJ • Additionally, 
the EOB Hamiltonian in [27| does not reduce to the Hamiltonian of a spinning test particle in the Kerr spacetime. 



A resolution to this issue was recently proposed in 31| , in which a canonical Hamiltonian was derived for a spinning 



test particle in a generic curved spacetime, to linear order in the particle spin. In turn, this canonical Hamiltonian 
has been used to construct an improved EOB Hamiltonian 32]. Additionally, the EOB model has recently been used 
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to model circular-equatorial extreme-mass-ratio inspirals (EMRIs) around spinning supermassive black holes, using 
fits of various post-Newtonian parameters to Teukolsky-based waveforms [331 ] - This is the same approach that was 
previously used with great success for non-spinning black hole systems, e.g., to derive fits for the final mass and spin 
of a BH after merger that are consistent with NR to about ~ 2% accuracy (34[. We used the non-spinning EOB 
model only in this work because that model is more mature. However, comparisons between an IMRI model based 
on the spinning EOB framework [33[ and the waveforms constructed using the transition model should be pursued in 
the future. 

The final ingredient, again common to both models, is the ringdown waveform. This originates from the distorted 
Kerr black hole formed after merger, and consists of a superposition of quasinormal modes (QNMs). Each mode 
has a complex frequency, whose real part is the oscillatory frequency, and whose imaginary part is the inverse of the 
damping time. Following Berti, et. al. [35{, and Buonanno, et. al. (36[, we developed a RD model that includes 
the fundamental mode (£ = 2, m = 2, n = 0) and two overtones (n — 1,2). QNMs always come in pairs which are 
called "twin modes" . In general, a mode with a given (£, m) will always contain a superposition of two different 
damped exponentials [35| , one with positive real part of the frequency, and the other with negative real part of the 
frequency and different damping time. Omitting these modes would have no serious consequences for non-rotating 
black holes, but it is conceptually inconsistent for rotating black holes. Furthermore, a single-mode expansion restricts 
attention to circularly polarized gravitational waves. For full generality, we use a general superposition of modes. The 
frequencies and damping times of the QNMs depend on the mass and spin after merg er of the newly formed Kerr 
BH and to estimate these we made use of fits for the final mass and spin found in [2l|, |37[ . The various QNMs are 
attached to the plunge and merger waveform at the effective light ring by imposing the continuity of the waveform 
and all necessary higher order time derivatives. 

Part of the motivation in developing a waveform model using the transition model of Ori & Thorne is to enable a 
more straightforward extension of this present waveform model for circular-equatorial IMRIs to a waveform model 
for generic (eccentric-inclined) IMRIs. The "numerical kludge" inspiral waveform [25[ and the transition model [38[ 
are both already available for generic orbits. The mass-ratio regime that we explore in these papers has not been as 
exhaustively explored as the comparable-mass-ratio and extreme-mass-ratio regimes. While the present paper was 
being prepared, results of a numerical simulation were published which shed light on the gravitational energy and 
momenta radiated by BH binaries with mass-ratios of 100:1 f39j . However, this was the first numerical simulation in 
the IMRI regime. The schemes described here combine elements of comparable mass and extreme-mass-ratio models 
in a logical way, but it cannot presently be tested against accurate simulations. As numerical and analytic models 
improve, these models should be refined. Nonetheless, the models should capture the main features of real IMRI 
waveforms and therefore be adequate for our purposes, namely to scope out the parameter estimation accuracies that 
might be achievable with ET. By using two alternative models we can build up confidence in the results. In [2(| we 
compared the two waveform models for non-spinning systems and demonstrated that the waveforms agreed well in 
both phase and amplitude, with differences at the level of a few percent. We show sample waveforms computed using 
the two models in Figure [2j The close agreement for non-spinning systems is clear from the figure in the lower panel. 
In this paper, we will now use these models to estimate the precision with which different third-generation detector 
networks will be able to estimate the parameters of the various binary systems listed above. 



In this section we will briefly review the aspects of signal analysis relevant to our calculations. A GW detector is a 
linear system whose input is the GWs we want to detect and whose output is a time series, which is a combination 
of the GW signals and instrumental noise. Any number of coplanar and colocated detectors have a response to GWs 
that can be written as a linear combination of the responses of two right-angle detectors offset by 45°. The output 
of each of these two equivalent two-arm Michelson detectors can be represented as 



where I, II denote the data streams from the two detectors. The detection problem is to distinguish h a (t) from n a (t). 
With the assumptions that (a) each Fourier component of the noise n a (f) has a Gaussian probability distribution; 
and (b) the noise is stationary, i.e., the different Fourier components are uncorrelated, the ensemble average of the 
Fourier components of the noise have the property 



IV. SIGNAL ANALYSIS 



s a (t) = h a {t) + n a {t), 



a = I, II, 



(1) 



(n a {f) Mf'T) = 2*(/ - f')S n (f)5 a p. 



(2) 



This relation defines the one-sided spectral density of the instrumental noise, S n (f). 
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FIG. 2: We show inspiral, merger and ringdown gravitational waveforms for compact objects of mass 10 Mq (top left panel and 
1.4M© (top right panel), inspiraling into a 500Mq BH with spin parameter q = 0.3. The bottom panel shows the gravitational 
waveform for a binary system consisting of a compact object of mass 10 Mq and a non-spinning (q — 0) IMBH of mass 5OOA/0, 
computed using both the EOB and transition models. In each case, the various extrinsic parameters were chosen randomly. 



We wish to estimate the parameters of a system from which we have detected gravitational waves. If the expected 
waveform h(t]9) depends on parameters 9 — #/v}, we want to construct the posterior probability of the 

parameters given the data and from that we will find the most probable value of the parameters of the source and 
their respective uncertainties. 

There is a natural inner product on the vector space of signals, which for any two signals p a (t),q a (t), takes the 
form 



(p|q)^2]T 



p*M)qM)+pM)q«(f) 

Sn(f) 



d/. 



The probability distribution for Gaussian noise, n(t), is given by 

( (n | n ) 
p(n ) = iVexp ' 



(3) 



(4) 



where N is a normalization factor. This relation gives the probability that the actual noise realization is n . If the 
source had parameters 9, the output would have the form s(t) = h(t; 9) + no(£). Using this relation and Equation (j4}, 
we obtain the likelihood of observing the output s(t) if the signal had parameters 9 



A(s\9\ = iVexp 



so -h(0)|s o -h(0) 



(5) 
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We can use Eq. ([5]) in Bayes' theorem, which states that the posterior probability is proportional to the product of the 
likelihood function and the prior probability. The prior describes the degree of belief, before a measurement is made, 
that a GW signal in the data would have parameters 9, and the posterior probability describes the corresponding 
degree of belief after the measurement. Denoting the prior by po(9), the posterior probability for the parameters 9, 
given the observed output s, is given by [40( 

p (h(0) | s) cx Po (0)exP ((V) |s) - \ (HO) | h(0))) . (6) 

As p (b-(0) \s^j is a distribution in 9 for a fixed output s, the term (s | s) /2 can be absorbed into the normalization 

factor. This probability distribution function describes the information we can extract from the data stream once the 
prior is specified. 

In the absence of an actual data set, we cannot compute this posterior probability distribution, but we can estimate 
how broad the posterior is likely to be and hence how well we might constrain the parameters. Assuming a flat 
prior (or a locally- flat prior which would be a reasonable assumption for a signal of high signal-to- noise ratio), the 
maximum posterior probability is the maximum-likelihood value, which minimizes (s — h | s — h). This is also the 
point in parameter space with the highest matched filtering signal-to-noise ratio (SNR) 



We can expand Eq. ([6]) about this peak, 9 by setting 9 l = 9 l + A9 l to obtain 



p(A0|*)« Ne^ T ^ A9 \ r 



ah 



dh 

w 



(8) 



Tij is the Fisher Information Matrix. This approximation is valid for large SNR. The correlation matrix of the errors 
A9 l is given by 

{A9 l A9 3 ) = (r- 1 )^. (9) 

We shall make use of Eq. (0) to estimate parameter errors for the various binaries described in Section III Bl For white 
noise, i.e., S n (f) = const., Parseval's theorem allows us to rewrite the inner product in the simple form |41| 



4s 



p a (t)q a (t)dt. (10) 



We shall make use of this alternative formulation of the inner product to simplify our analysis. If we define the 
"noise-weighted" waveform 

Kit) - -tMU /(*) = I§, (11) 
^ (/(*)) 



we can then rewrite the Fisher matrix approximately as |4l| 

l-T 



r Qb = 2^/ d a h a (t)d b h a (t)dt , where d a = d/d6 a . 

r. JO 



(12) 



As discussed in [20} , the binaries we shall consider have relatively low SNR and it is known that in such situations 
the Fisher Matrix may overestimate measurement accuracies. Vallisneri [42| gave a mismatch criterion to determine 
whether Fisher Matrix results were reliable in a particular context. The idea is to compute the ratio r(9, A) of the 
linearized signal amplitude (LSA) likelihood to the exact likelihood. This ratio r is given by [42] 

|logr(0,A)| = (AOjhj ~ Ah{9),A9 k h k - Ah(9)) , (13) 

where A represents the signal strength, Ah(9) = h(9) — h(9) , with 9 = 9 + A9, and 9 the observed location of maximum 
LSA likelihood for a given experiment. We use the consistency criterion, Eq. (fl"3j) . to explore the \-a likelihood surface 
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predicted by the Fisher Matrix to verify that the mismatch between the LSA and the exact likelihoods is small. Ratios 
below a fiducial value, say | logr(0, A)| ~ 0.1, are considered acceptable [42j |. 

As an example of typical results obtained from this analysis, for (1OM , 100M Q ) binaries with q — 0.3, we found that 
systems near the lower quartile of the distribution had | \ogr(8, A)| ~ 0.2. From this threshold onwards, the ratio r 
decreases gradually so that at the upper quartile of the distribution we obtain | \ogr(9, A)| ~ 0.04. This indicates that 
the results we get from the Fisher Matrix should be a reasonable estimate of the measurement precisions achievable 
for IMRIs using ET. Our results may be somewhat optimistic and at some point should be verified using Monte Carlo 
simulations to recover the full posterior probability distributions. Such an exercise is beyond the scope of this present 
work. Nonetheless, the results should be fairly accurate, and we will see in Section [V] that, where comparisons can be 
made, our results are in good accord with previous results that have been derived for other types of binary and using 
different waveform models. 

The Fisher matrix for a network of detectors is given by the sum of the individual Fisher Matrices for each detector. 
When modelling the response of the various detectors, it is necessary to account for their relative positions on the 
surface of the Earth, as the corresponding time delays are what allow source triangulation. In computing the waveforms 
used in the Fisher Matrix, we found it convenient to use two different timesteps in order to separately resolve the 
(slow) inspiral and (fast) merger/ringdown phases. We checked that varying the choice of timesteps and the location 
of the transition between the two timesteps did not significantly affect the results for any of the systems considered. 
We will discuss the convergence of the Fisher Matrices further in Section fVl 

V. PARAMETER ESTIMATION ERROR RESULTS 

The parameter space of the signals we shall consider is 10 dimensional. Four of these are intrinsic parameters, 
namely In m, InM, q, to. The other six are extrinsic parameters. We summarize the physical meaning of all the 
parameters in Table HI 

To explore parameter estimation errors using the inverse Fisher Matrix we fix the values of the intrinsic parameters 
of the source, and carry out a Monte Carlo simulation over possible values for the extrinsic parameters. The parameter 
errors scale with the SNR of the source as 1/SNR, so we can quote results at a fixed SNR and these may be easily 
extrapolated to other SNRs. We chose a reference SNR of 30, as this would represent a very robust detection for 
which the Fisher Matrix prediction is likely to be a good estimate. To obtain this reference result, we first compute 
the Fisher Matrix for a source at a fixed distance, D = 6.6348Gpc, and obtain the SNR at that distance from the 
expression 




We then multiply the errors estimated from the Fisher Matrix by (SNR/30) to normalise to the reference SNR of 30. 
We assume that the observation starts when the GWs from the inspiral reach a frequency of 5Hz and finishes when 
the RD waveform is no longer contributing to the SNR. 



lnm 


mass of compact object (CO) 


InM 


mass of IMBH 


q 


magnitude of (specific) spin angular momentum of IMBH 


to 


time at which orbital frequency sweeps through a reference value 


00 


initial phase of CO orbit 


e s 


source sky colatitude in an ecliptic-based system 


4>s 


source sky azimuth in an ecliptic-based system 


6k 


direction of IMBH spin (colatitude) 




direction of IMBH spin (azimuth) 


InD 


distance to source 



TABLE I: This table shows the physical meaning of the parameters used in our model. The angles (8s, <fis) an d (@Ki 4>k) are defined in 
a fixed ecliptic— based coordinate system. 
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A. Dependence of parameter estimation errors on system parameters 

We consider the twelve different binary systems described in Section III Bl For the sources with q = 0.3, 0.9 
we present results computed from the transition model for the waveform. For the non-spinning systems, q = 0, 
we present results from both the transition model and the EOB model. The Tables HlHVl list the mean, standard 
deviation, median and lower and upper quartiles of the distribution of the Fisher Matrix errors computed in the 
Monte Carlo simulation. There is one table for each of the four mass combinations we consider and all results are 
computed at the fixed reference SNR of 30. These results assume that the detection is made using the optimistic 3 
ET network configuration, and hence should be considered as upper bounds on the accuracy with which we may be 
able to measure the various parameters. 

In Figure [3] we show histograms of the error estimates for several parameters found in the Monte Carlo simulation 
for the m = WMq, M — lOOAf©, q — 0.9 system. This figure, and Table ITT), indicate that a network of 3 ETs might be 
able to determine the location of the source in the sky to an accuracy of ~ 10~ 3 steradians, i.e., ~ 4 square degrees. 
This estimate is the statistical mean of the Monte Carlo distribution at SNR of 30 assuming the optimistic detector 
network of three ETs in the geographic locations of Virgo, LIGO Livingston and Perth, Australia (AIGO). 

Recent work has shown that for the existing LIGO- Virgo detector network, assuming a uniform distribution of 
sources across the sky, at a network SNR of around 15, 50% of inspiral sources should be located within 23 sq-degs 
(best case) at the 95% confidence level. For burst sources, without any knowledge of the waveform, 50% of the sources 
could be localized within 50 sq-degs (worst case) at an SNR of 10, but this can be reduced to 8 sq-degs if predicted 
waveforms are available These estimates are similar for the initial or advanced detectors (43|], at a fixed SNR. 
The inclusion of an additional detector in the Southern Hemisphere, such as AIGO, can further improve these values 
as it contributes a longer baseline, additional energy flux, extended signal space, and breaks the plane-degeneracy 
formed by the three detector network of LIGO Livingston, LIGO Hanford and Virgo [43^ . 

Another recent independent study [44[ presented results for the accuracy with which sources could be localized 
with a network of GW detectors, using only timing information in the various detectors. Fairhurst [44| found that 
increasing the number of detectors at different sites increases both the absolute number of observable sources, and 
greatly increases the fraction of sources that can be well localized. For instance, at a fixed SNR of 8, and using the 
advanced detector network comprising the two LIGO detectors at Hanford (HH) and LIGO Livingston (L), no sources 
can be localized within 20 sq-degs. Adding another detector at a different site, e.g., a LIGO detector in Australia 
(A), or Advanced Virgo (V), or the Japanese LCGT detector (J), the networks HLA, HHLV, HHJL, can localize up to 
50% of the signals within 20 sq-degs, and the loudest signals within 5 sq-degs. In all these results, it was found that 
for the networks involving an Australian detector, the peak of the localization distributions occurs between 5 and 10 
sq-degs. Our results imply that we can determine the location of the source in the sky to an accuracy of ~ 12 sq-degs 
when we renormalise to a network SNR of 8. This estimate is the statistical median value obtained from the Monte 
Carlo simulations. Our results are therefore in good accord with existing estimates of the angular resolution that 
could be achieved by the advanced detector network. Furthermore, Table |TT] shows that for non-spinning systems, we 
should resolve the plunge time to to within ~ 30 ms when normalised to an SNR of 10. This estimate is also in good 
accord with the results obtained independently by Ajith et al. [45| and Luna et al. [46| . 

The Monte Carlo simulations also show that we expect to determine the total mass of m = WMq , M — 100M Q 
binaries to an accuracy of ~ 0.1% for non-spinning systems when we rcnomalise the results to an SNR of 8. This 
estimate is in good accord with previous estimates for related systems. In [47} it was estimated that using the same 
3ET detector network that we have considered in this paper, it should be possible to determine the total mass of a 
non=spinning m = 23M Q , M = 100M© binary (mass ratio rj = 0.16) to an accuracy of ~ 0.1% at a network SNR 
of 8. This estimate was obtained using the phenomenological waveform model described in [48j . which includes the 
inspiral, merger and ringdown for non-spinning, comparable-mass binaries in a consistent way. 

Figure 2] shows the accuracy with which we expect to determine the masses of the small CO and the IMBH for the 
spinning systems, as well as the IMBH's spin parameter q. This figure shows that for binaries that stay the longest in 
band, namely the m = 1.4M© , M = WQMq, q — 0.3 system ( ~ 291 seconds in band), we can determine the masses 
of the CO and the IMBH to an accuracy of 0.01%, 0.005%, respectively, at the reference SNR of 30. The accuracy 
gradually decreases for the shorter-lived events. For instance, for binaries with m = 10M Q , M = 500M Q , q — 0.3, 
which stay in band about ~ 8.3 seconds, we expect to measure the CO's and IMBH's masses to an accuracy of only 
- 0.15%, - 0.05% respectively, at the same SNR of 30. 

We expect to determine the spin parameter more accurately for rapidly rotating binaries. For instance, the estimated 
accuracy with which we can measure the IMBH's spin for binaries with m — 10M Q , M = WOMq, q — 0.9 is ~ 0.05%. 
This is about a factor of 2 better than for binaries with the same component masses but with spin q — 0.3. We would 
expect this trend since, for more rapidly spinning systems, the small CO comes much closer to the outer horizon of 
the Kerr IMBH before merging, and this is the regime where the CO can more strongly feel the effects of the IMBH's 
spin. We note that the percentage error in the spin is quite similar in the two cases. This is consistent with this 
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FIG. 3: Distribution of parameter measurement error estimates computed in the Monte Carlo simulations. We show results 
for the system with m — IOA/q, M = IOOM0, q — 0.9. The panels show the error distributions in the following order, top row, 
from left to right, A(lnm), A(lnM); middle row, Aq, Ato; bottom row, AOs, A(j>s. 
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FIG. 4: Left panel: expected fractional measurement errors, A, in the mass of the CO and the mass of the IMBH for a fixed 
value of the IMBH spin parameter q = 0.3. There is a pair of candlesticks for each of the four binary systems. For each binary, 
the candlestick representing the CO's mass error estimate is to the left of that for the IMBH's mass error estimate. The error 
bars indicate the lower and upper quartiles found in the Monte Carlo simulation. Right panel: expected measurement errors in 
the IMBH's spin parameter for two different values for the spin parameter q. For each system, the candlestick corresponding to 
a central IMBH with q = 0.3 is to the right of that associated with the error estimate for a system with q = 0.9. The continuous, 
long-dashed, dashed and dotted lines represent systems Bl: m = 10M Q ,M = 100M©, B2: m = 1.4M Q ,M = 100M Q , B3: 
m = 10M Q , M = 500M o , and B4: m = 1.4M© , M = 5OOM , respectively. 



understanding — the spin is measured more precisely when it is larger because it has a correspondingly greater effect 
on the system. For less massive COs, the inspiral proceeds more slowly and the early inspiral has a correspondingly 
greater importance in the estimation of the system parameters. For such systems, we might expect the precision of 
spin determination to be les sensitive to the system spin, and this is borne out by our results — the spin precision is 
comparable for q = 0.9 and q = 0.3 in the 1.4M© + 100M Q and 1.4M + 5OOM systems. 

In general, we expect that the precision with which we can measure the parameters of a binary depends on the 
number of GW cycles that are observed. Therefore, we might expect the precision to depend on the source masses 
as, from best to worst, 1.4M + 1OOM , 1.4M Q + 500M Q , 1OM + 100M o , 1OM + 500M Q . Tables HUd and 
Figure U confirms this expectation. These Tables also show that the precision improves for more rapidly spinning 
IMBHs. This is because, for greater spins, the inspiral phase evolution lasts longer and so there are more cycles of 
information in the waveform. The SNR of the spinning systems will also be larger at a given distance, as compared 
with slowly spinning IMBHs, so these systems can be seen further away. This does not affect the results quoted here, 
which are normalised to fixed SNR=30. Tables iLTHVl also confirm that our two independent waveform models make 
predictions for the parameter estimation errors for non-spinning systems that are consistent to better than ten per 
cent. 

Finally, we note that the consistency of the distributions of Figure |3l in the sense that they are smooth with 
few outliers, is an indication that the results are convergent. It is well known that Fisher Matrices encountered in 
parameter estimation calculations for GW sources can have very large condition numbers. However, the results we 
present here were obtained from Fisher Matrices that exhibited convergence over at least two orders of magnitude 
in the offsets used to compute the numerical waveform derivatives. The inverse matrices were computed using an 
LU decomposition, and we verified that the inverse Fisher matrices were also convergent to <J 10% over an order of 
magnitude in the numerical offsets. We also found that the offsets required for the various network configurations 
were consistent, as we would expect, since the convergence of the FM should depend on the intrinsic waveform, rather 
than the choice of network. 



B. Dependence of parameter estimation errors on network configuration 



We shall now consider the more modest network configurations C1-C4 described earlier, which consist of combi- 
nations of ET's and right-angle detectors at 2 or 3 sites chosen from Virgo, LIGO Livingston and Perth (Australia). 
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These results are shown in Tables IVTrllXI 1 for the four different mass combinations, but with the central IMBH spin 
fixed at q = 0.3. As expected, we find that a single ET is sufficient for accurate intrinsic parameter determination. 
This is because the signature of the intrinsic parameters is encoded in the phase evolution which can be accurately 
measured by a single instrument. After normalising the parameter errors to a fixed network SNR of 30, Figure [5] 
shows clearly that the accuracy with which we can measure the intrinsic parameters of a binary is not dramatically 
changed by changes in the network configuration. 



0.001 - 




A 



0.0001 1 1 

FIG. 5: Arranged from left to right in the panel, we show the expected measurement errors in the CO's mass, IMBH mass and 
IMBH spin parameter, as a function of the network configuration for binary systems with masses IOMq + WOMq and IMBH 
spin q — 0.3. In each group the candlesticks are arranged in order C1-C5 from left to right. As before, the CO and IMBH 
mass errors are quoted as fractional errors, i.e., these are the errors in ln(m), ln(Af). 

However, it is not possible to constrain the location of the source in the sky nor the luminosity distance using a 
single detector. This is also expected, as the determination of the extrinsic parameters comes primarily from the 
time-delays between the arrival of the signal in different detectors. It appears that the inclusion of an additional 
right-angle detector (configuration C2) is enough to constrain the source sky position and luminosity distance to 
moderate precision ~ 30 sq-degs, ~ 15% for a 10M© BH IMRI into a 100M Q IMBH with spin q = 0.3 at SNR of 30. 
Replacing this right-angle detector by another ET at the same location (configuration C3) enhances the determination 
of these parameters to ~ 12 sq-degs, ~ 10%, respectively. If we consider a configuration consisting of 1 ET, plus two 
right-angle detectors (configuration C4) these estimates are further improved to ~ 11 sq-degs, ~ 10%, respectively. 
These two estimates are rather close to the precision we expect to obtain using a 3 ET network, namely, ~ 8 sq-degs, 
~ 10%, respectively. 

Moreover, these results also show that the distributions become narrower for the more complex configurations, i.e., 
there is less variation in the precisions with which we can determine the extrinsic parameters as we randomise over 
the source position and orientation. This is also expected, as a more complex network should provide more complete 
sky coverage. We can compare our results for the luminosity distance in the non-spinning limit with those quoted 
in Table I of (47[, which were given at a network SNR of 8. Renormalising our results for the non-spinning systems 
to an SNR of 8, we estimate that network C2 would be able to measure the luminosity distance to an accuracy of 
~ 40%. For configuration C3, we would be able to determine the luminosity distance to an accuracy of 30% and for 
C4, this improves marginally to ~ 28%. This estimate is roughly the same for the 3ET detector network (C5). These 
estimates are in good accord with Table I of |47|, which were computed for a ~ 20M Q + 100M Q binary system using 
a comparable- mass inspiral-merger-ringdown waveform model. This comparison provides further support that even if 
the IMRI models used in this paper are not accurate enough to be used for detection of IMRI systems, they do still 
give reasonable results and can be used for parameter estimation studies to illustrate the potential scientific impact 
of the ET. 

In summary, a single ET will be enough to measure the intrinsic parameters of a system but will not be sufficient 
to accurately reconstruct the extrinsic parameters. Extrinsic parameter determination is possible with the addition of 



1 Note that only Table IVT1 is included in the main text. Tables IVimiXI may be found in Appendix A. 
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one right-angle detector, but a 2 ET configuration will improve the accuracy with which we can measure the source's 
sky position and luminosity distance by roughly a factor of ~ 2. A more optimistic configuration consisting of a single 
ET and two right-angle detectors upgraded to ET's sensitivity will generate results that are competitive with the 
highly optimistic three ET network. This configuration (C3) might be more realistic, since it will have lower overall 
costs, although it is not simply an upgrade of existing sites as we are assuming 10km scale right-angle interferometers. 

Tables iHllIXI show that GW observations will produce valuable astrophysical results. A network of 2 ETs operating 
in coincidence will be able to measure accurately the masses of two merging black holes, and less accurately the sky 
position and luminosity distance at which the merger is taking place. The results in paper I [20j suggested that a 
a few tens to a few hundreds of IMRI events could be detected per year by an ET network, provided that IMBHs 
are relatively abundant in globular clusters, and the various mechanisms leading to CO capture — secular Kozai 
resonance, binary exchange processes, gravitational radiation and three and four-body interactions — are efficient at 
starting IMRIs. The ET network would be able to achieve good parameter determinations for all the systems that 
are detected, which will have important scientific implications. 

One important question will be whether the IMBHs detected have been formed in clusters, or are primordial, 
population III, black holes formed in the early Universe. The precision achievable with network C4, of ~ 11 sq-degs 
in sky position and ~ 10% in luminosity distance will not be enough to determine the host galaxy of an IMRI uniquely, 
let alone if it is occurring in a cluster within that galaxy. If the the ET low-frequency cut-off does extend down to 
1Hz, these precisions improve to ~ 1 sq-degs, and ~ 7%, but even that is unlikely to be enough. However, primordial 
IMBHs are likely to be hosted in dwarf galaxies today, and the stellar density is sufficiently low that is very unlikely 
that IMRIs could occur there [ll|. Thus, it would be a reasonable assumption that any IMRIs detected would be from 
cluster IMBHs and not population III black holes. The measured masses and IMBH spins thus provide constraints 
on the cluster IMBH population, and their evolution. 

IMRI systems could be detected up to redshifts as high as z ~ 6, if the system has redshifted masses [10 + 1OO]M 
and IMBH spin q = 0.9, assuming a 3ET network with a frequency cut-off of 1Hz 20]. At such redshifts, the existence 
of light black holes, no matter how they formed, provides some constraints on the hierarchical assembly of structure. 
If IMBHs exist at z ~ 6, they will inevitably grow and may be progenitors of some fraction of today's supermassive 
black holes. IMRI detections with ET could therefore be important for cosmology, which relies on the ability to 
confidently say that an event is at high redshift, which we have demonstrated is possible. 

A redshift z ~ 6 is probably not high enough to provide strong constraints on structure growth, but ET might also 
detect IMBHs when they merge with comparable mass IMBHs, and these can be seen out to redshifts z J> 10 fill. [47l| . 
Again, the IMBHs involved in such mergers could form either in clusters or in the early Universe. As many as a few 
tens of mergers involving primordial IMBHs could be detected by ET [47| in the light-seed scenario. The number 
of comparable mass mergers of IMBHs formed in clusters could be much higher, as many as several thousand fTTj | . 
but this is very dependent on the efficiency of IMBH formation in the cluster environment. As mentioned above, 
the IMRIs detected by ET are most likely to be occurring in clusters and thus provide a direct constraint on the 
existence, number density and properties of IMBHs in clusters. This information can be used to estimate how many 
of the observed IMBH-IMBH mergers might be coming from the cluster channel and hence what fraction of events 
might be primordial. IMRI observations will therefore also have indirect applications to our understanding of structure 
formation. 

The Monte Carlo results we present in this paper are the first results for IMRI sources detectable by ET to appear 
in the literature. The results must be taken with some caution, because of the approximations in the waveform models 
that have been made, as discussed earlier. However, these results should be a reasonable guide to the likely order-of- 
magnitude of the parameter measurement errors that will be achievable. The precise astrophysical implications will 
depend primarily on the number of events that are seen, which is very difficult to predict given current uncertainties 
in the astrophysics of IMBHs, and indeed their very existence. 

VI. CONCLUSIONS 

In this paper, we have used the gravitational waveform models for IMRIs developed in paper I of this series [2(| to 
estimate the precision with which the Einstein Telescope will be able to determine the parameters of circular-equatorial 
IMRIs. We have presented results for a set of twelve "typical" systems, comprising four different combinations of 
component masses — 1.4M© + 100M o , 1.4M + 5OOM , 1OM + 100M Q , and 1OM + 500M o — and three different 
IMBH spins — 9 = 0, 0.3, 0.9. For the non-spinning systems, we have compared the results from the transition model 
and the EOB waveform model and found that these models make predictions that are consistent to better than ten 
percent. This final check provides confidence in these results. 

We have also explored how the accuracy of parameter determination depends on the configuration of the detector 
network using the "ET-B" noise curve and assuming a cut-off frequency of 5Hz. We have shown that a single ET 
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TABLE II: Summary of Monte Carlo results for parameter estimation errors. We show the mean, standard deviation, median 
and quartiles of the distribution of the error in each parameter. Results are given for a m = IQMq CO inspiralling into a 
M = IOOMq IMBH for various choices of the IMBH spin, q, computed using the transition model waveform. We also show 
results for q = computed using the EOB waveform model. 

is sufficient to accurately determine the intrinsic parameters of these systems. However, a network of detectors is 
required to obtain accurate estimates of the extrinsic parameters. A network of 2 ETs should be sufficient to measure 
the source's sky position and luminosity distance to accuracies of ~ 12 sq-degs and ~ 10%, respectively for a BH 
IMRI into a IOOM0 IMBH with spin q = 0.3. A more sophisticated network comprising 1 ET and two right-angle 
detectors would have comparable precisions. The results for the latter networks are comparable to the precisions that 
could be achieved with a 3 ET network, which are ~ 8 sq-degs, 10%, respectively (at a source SNR of 30). Any one of 
these ET networks would simultaneously be able to constrain the BH and IMBH masses and the IMBH spin magnitude 
to fractional accuracies of ~ 10~ 3 , 10~ 35 and 10 -3 , respectively. The amount of variation in the parameter precision 
over random choices of the source location and orientation also decreases for more complex network configurations. 
A 3 ET network is a highly optimistic assumption about a future third-generation GW detector network, but our 
results indicate that a more modest network comprising one ET and right-angle interferometers in LIGO Livingston 
and Perth can recover parameters to almost the same precision. This network would have lower associated costs 
and might therefore be more feasible. Our results should be regarded as conservative in the sense that using a lower 
low- frequency cut-off, or assuming a different ET design, e.g., the xylophone configuration [i"4| , could improve the 
various results quoted in these papers. At present we have made use of the "ET-B" noise curve (23[, but a study of 
the potential applications of ET using both lower cut-off frequencies and more optimistic designs should be carried 
out in the future. 

If IMRIs are detected, there will be several important scientific payoffs. The detection of an IMBH as a GW source 
would provide the first robust proof of the existence of black holes in the intermediate mass range. As discussed 
in the previous section, IMRIs are most likely to occur in globular clusters, thus IMRI observations will provide 
information about whether IMBHs form at all in these stellar environments, and whether they remain in the clusters 
rather than being ejected. The masses, spins and abundance of IMBHs detected as IMRIs will provide constraints on 
the formation efficiency and evolution of globular cluster systems. In conjunction with ET observations of comparable 
mass IMBH-IMBH mergers, IMRIs could shed light on the hierarchical growth of structure, in particular to distinguish 
between light and heavy seed models, as described in more detail at the end of Section IVl 

A further potentially exciting payoff of ET IMRI observations will be to test whether the central object in an 
IMRI is described by the Kerr metric of general relativity. During the inspiral, the CO traces out the geometry of the 
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Statistics of distribution for error AX in parameter X = 
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TABLE III: As Table Ull but now for binary systems with a CO of mass m = 1.4Mq. 



spacetime of the central object. Hence, the GWs emitted during the inspiral encode a map of the spacetime. In [49l |50|| 
it was shown that any stationary, axisymmetric, vacuum spacetime in relativity can be decomposed into mass, Mi, 
and current, St, multipole moments. Ryan [5l[ went on to show that, for nearly equatorial, nearly circular inspirals, 
these multipole moments are redundantly encoded in gravitational wave observables, namely the periapsis precession 
frequency, the orbital plane precession frequency and the gravitational wave energy spectrum. The multipole moments 
of a Kerr black hole are completely determined by its mass, M, and spin, Si, through the relation [50j 

M e + iSi = M (iaf , (15) 

where a — S\/M is the reduced spin of the black hole. Extracting three moments of the spacetime from GW emission, 
and finding them inconsistent with (fT"5|) . would suffice to demonstrate that the central object is not a Kerr black hole. 
Tests of this nature have been considered for observations of EMRI systems with LISA. Using a kludge model that 
included a non-Kerr value for the quadrupole moment of the central black hole, Barack and Cutler |52i | showed that 
LISA could measure the quadrupole moment, Q = —S 2 /M, of the central black hole to an accuracy AQ/M 3 ~ 10~ 3 , 
while simultaneously measuring the mass and spin to an accuracy ~ 10~ 4 . In the context of IMRIs detected by 
Advanced LIGO, it has been shown that it is possible to measure an 0(1) fractional deviation in the mass quadrupole 
moment for typical systems 53]. ET could improve this precision by more than an order of magnitude, since ET's 
enhanced sensitivity will allow a significantly larger number of gravitational-wave cycles to be observed than will be 
possible with Advanced LIGO. For instance, for a IMq + 100M© system Advanced LIGO could measure ~ 500 cycles 
until plunge. In contrast, ET with a low-frequency cut-off at 1Hz will observed up to ~ 25000 cycles [H}. Further 
work is required to properly quantify the ability of ET to carry out these sorts of tests on the nature of IMBHs, and 
this should be pursued in the future. 

Our results are a first attempt to explore the precision of IMRI parameter estimation that will be achievable with the 
Einstein Telescope. Our results have been derived using a particular waveform model, which reflects the best of what 
is currently available, and combines results from both comparable mass binaries and extreme-mass-ratio inspirals. 
These waveforms are unlikely to be accurate enough to be used in a search to recover source parameters, but they 
should capture most of the main features of true IMRI waveforms and therefore provide a decent estimate of the 
level of precision that could be achieved by GW measurements. The waveform models could be improved in various 
ways, by including conservative corrections, by generalizing the waveform models to consider both eccentric orbits 
and orbits inclined to the equatorial plane, and by including the leading order effects of the spin of the smaller object. 
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Statistics of distribution for error AX in parameter X = 
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TABLE IV: As Table HH but now for binary systems with a central IMBH of mass M = 500M©. 



The EOB model can also now be extended to spinning systems, for circular-equatorial inspirals at least |33(, which 
will provide further important consistency checks for the transition model. These will be important improvements to 
consider in the future in order to confirm the present results, and to extend the calculations to generic IMRI systems. 
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Statistics of distribution for error AX in parameter X = 



Model 


ln(m) 


ln(M) 


q 


to 


</>o 


e s 


<t>s 


9 K 


4>K 


ln(D) 




Mean 


1.23e-3 


3.45e-4 


1.15e-3 


2.04e-2 


1.62e-l 


2.2995 


2.7421 


3.7537 


4.0738 


4.2386 




St. Dev. 


4.19e-4 


1.99e-4 


4.23e-4 


9.71e-3 


5.90e-2 


3.4104 


4.0524 


3.9352 


4.0112 


4.1351 


CI 


L.Qt. 


1.02e-3 


2.01e-4 


8.12e-4 


1.41e-2 


1.23e-l 


0.0866 


0.1018 


1.0189 


1.1825 


1.2491 




Med. 


1.24e-3 


3.07e-4 


1.15e-3 


1.99e-2 


1.62e-l 


0.6820 


0.8174 


2.2722 


2.6118 


2.6707 




U. Qt. 


1.43e-3 


4.55e-4 


1.55e-3 


2.81e-2 


1.84e-l 


2.9151 


3.3087 


4.9711 


5.3026 


5.7026 




Mean 


l.lle-3 


3.13e-4 


9.12e-4 


1.69e-2 


1.34e-l 


6.19e-2 


7.61e-2 


1.50e-l 


1.88e-l 


1.76e-l 




St. Dev. 


4.84e-4 


1.84e-4 


3.57e-4 


8.43e-3 


4.35e-2 


3.77e-2 


4.62e-2 


1.28e-l 


1.76e-l 


1.54e-l 


C2 


L.Qt. 


9.43e-4 


1.81e-4 


7.41e-4 


1.23e-2 


1.08e-l 


3.31e-2 


3.76e-2 


6.60e-2 


6.34e-2 


9.05e-2 




Med. 


1.15e-3 


2.72e-4 


9.77e-4 


1.82e-2 


1.45e-l 


5.18e-2 


6.22e-2 


1.01e-l 


1.28e-l 


1.47e-l 




U. Qt. 


1.33e-3 


4.19e-4 


1.26e-3 


2.63e-2 


1.75e-l 


7.79e-2 


1.07e-l 


2.03e-l 


2.48e-l 


2.74e-l 




Mean 


1.07e-3 


3.04e-4 


8.87e-4 


1.66e-2 


1.29e-l 


4.56e-2 


5.81e-2 


8.47e-2 


1.14e-l 


l.lle-1 




St. Dev. 


3.73e-4 


1.75e-4 


3.59e-4 


8.01e-3 


4.12e-2 


2.15e-2 


3.27e-2 


6.78e-2 


9.85e-2 


8.35e-2 


C3 


L.Qt. 


9.09e-4 


1.78e-4 


6.91e-4 


1.12e-2 


1.08e-l 


2.89e-2 


3.24e-2 


4.07e-2 


4.92e-2 


5.81e-2 




Med. 


l.lle-3 


2.64e-4 


9.12e-4 


1.77e-2 


1.45e-l 


4.29e-2 


4.92e-2 


6.07e-2 


7.83e-2 


8.73e-2 




U. Qt. 


1.29e-3 


3.98e-4 


1.26e-3 


2.63e-2 


1.75e-l 


5.95e-2 


7.20e-2 


9.78e-2 


1.49e-l 


1.45e-l 




Mean 


1.05e-3 


3.01e-4 


8.79e-4 


1.58e-2 


1.18e-l 


4.47e-2 


5.61e-2 


8.24e-2 


1.14e-l 


1.12e-l 




St. Dev. 


4.27e-4 


1.78e-4 


3.56e-4 


8.14e-3 


4.38e-2 


2.09e-2 


2.87e-2 


6.38e-2 


9.81e-2 


8.07e-2 


C4 


L.Qt. 


8.38e-4 


1.75e-4 


6.31e-4 


1.12e-2 


1.04e-l 


3.21e-2 


3.54e-2 


4.46e-2 


4.82e-2 


5.38e-2 




Med. 


1.12e-3 


2.69e-4 


8.71e-4 


1.77e-2 


1.45e-l 


4.06e-2 


5.02e-2 


6.91e-2 


7.71e-2 


8.52e-2 




U. Qt. 


1.27e-3 


4.02e-4 


1.20e-3 


2.39e-2 


1.69e-l 


5.45e-2 


6.91e-2 


9.54e-2 


1.33e-l 


1.40e-l 




Mean 


1.09e-3 


2.96e-4 


8.76e-4 


1.51e-2 


1.13e-l 


4.28e-2 


5.42e-2 


7.80e-2 


9.77e-2 


1.05e-l 




St. Dev. 


3.68e-4 


1.22e-4 


3.13e-4 


8.39e-3 


3.55e-2 


1.54e-2 


2.48e-2 


5.17e-2 


4.35e-2 


7.14e-2 


C5 


L.Qt. 


8.93e-4 


1.73e-4 


6.31e-4 


1.12e-2 


8.39e-2 


2.95e-2 


3.54e-2 


3.93e-2 


4.64e-2 


5.53e-2 




Med. 


l.lle-3 


2.64e-4 


8.51e-4 


1.66e-2 


1.29e-l 


4.16e-2 


5.04e-2 


5.73e-2 


7.45e-2 


8.64e-2 




U. Qt. 


1.22e-3 


3.91e-4 


1.01e-3 


2.34e-2 


1.59e-l 


5.22e-2 


6.96e-2 


8.64e-2 


1.24e-l 


1.39e-l 



TABLE VI: As Table ITU but for binary systems with a central IMBH of spin parameter q = 0.3 and assuming four alternative 
configurations for the detector network, C1-C4, as described in Section [TTJ Configuration C5 is the network of three ETs which 
has been used for all results elsewhere in this paper. 
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Appendix A: Dependence of parameter estimation errors on network configuration 

This Appendix contains Tables IVIIHIX1 These Tables show how the expected parameter measurement errors 
depend on the network configuration. We explore how parameter estimation accuracies are modified for five network 
configurations, C1-C5. These configurations are, CI: one ET at the geographic location of Virgo; C2: as configuration 
CI plus a right-angle detector at the location of LIGO Livingston; C3: as configuration CI plus another ET at the 
location of LIGO Livingston; and C4: as configuration C2 plus another right-angle detector in Perth. We denote the 
reference 3-ET network as configuration C5. 



Model 


Statistics of distribution for error AX in parameter X — 


ln(?77.) 


ln(M) 


n 

y 








<i>s 


0k 




In CD") 


CI 


Mpan 


1.26e-4 


7.76e-5 


1.96e-4 


2.63e-2 


2.01e-l 


3.7418 


4.1009 


6 6255 


7.4131 


9.2443 


St. Dev. 


6.54e-5 


4.57e-5 


9.27e-5 


1.30e-2 


1.09e-l 


4.0925 


4.7919 


4.5011 


5.2979 


6.6406 


L.Qt. 


6.73e-5 


4.47e-5 


1.84e-4 


1.73e-2 


1.29e-l 


0.5762 


0.4917 


2.8422 


3.5419 


5.2586 


Med. 


9.93e-5 


7.08e-5 


2.19e-4 


2.51e-2 


1.85e-l 


2.0259 


2.0586 


5.8983 


6.9206 


8.5296 


U. Qt. 


1.66e-4 


1.17e-4 


2.36e-4 


3.54e-2 


2.56e-l 


6.1471 


6.1689 


8.9838 


11.4924 


12.6824 


C2 


Mpan 


1.09e-4 


6.61e-5 


1.84e-4 


1.77e-2 


1.64e-l 


5 85e-2 


7.58e-2 


1.70e-l 


2.34e-l 


1.98e-l 


St. Dev. 


5.84e-5 


3.73e-5 


9.13e-5 


1.02e-2 


8.48e-2 


3 08p-2 


4.71e-2 


1.40e-l 


2.78e-l 


1.53e-l 


L.Qt. 


6.17e-5 


4.07e-5 


1.44e-4 


1.34e-2 


1.14e-l 


3.30e-2 


3.56e-2 


6.64e-2 


7.24e-2 


8.31e-2 


Med. 


8.81e-5 


6.17e-5 


1.90e-4 


1.65e-2 


1.50e-l 


5.22e-2 


6.02e-2 


1.28e-l 


1.81e-l 


1.45e-l 


U. Qt. 


1.46e-4 


1.02e-4 


2.19e-4 


2.33e-2 


2.17e-l 


7.73e-2 


1.04e-l 


2.36e-l 


5.01e-l 


2.50e-l 


C3 


Mean 


9 82e-5 


6.31e-5 


1.81e-4 


1.65e-2 


1.50e-l 


4.36e-2 


4.95e-2 


8.71e-2 


1.44e-l 


1.05e-l 


St. Dev. 


5.28e-5 


3.48e-5 


9.09e-5 


9.11e-3 


7.71e-2 


1.55e-2 


2.15e-2 


6.34e-2 


1.39e-l 


7.11e-2 


L.Qt. 


5.62e-5 


3.80e-5 


1.29e-4 


1.28e-2 


1.17e-l 


2.95e-2 


3.06e-2 


4.36e-2 


5.37e-2 


5.46e-2 


Med. 


8.11e-5 


5.62e-5 


1.86e-4 


1.54e-2 


1.47e-l 


4.57e-2 


4.49e-2 


6.30e-2 


1.28e-l 


7.94e-2 


U. Qt. 


1.31e-4 


9.77e-5 


2.14e-4 


2.13e-2 


1.95e-l 


5.72e-2 


6.54 e-2 


1.14e-l 


2.69e-l 


1.31e-l 


C4 


Mean 


9.79e-5 


6.03e-5 


1.81e-4 


1.65e-2 


1.50e-l 


3.97e-2 


4.67e-2 


8.74e-2 


1.41e-l 


1.02e-l 


St. Dev. 


4.78e-5 


3.37e-5 


8.81e-5 


8.81e-3 


5.84e-2 


1.18e-2 


1.72e-2 


6.33e-2 


1.20e-l 


6.84e-2 


L.Qt. 


5.66e-5 


3.74e-5 


1.29e-4 


1.31e-2 


1.17e-l 


3.03e-2 


3.17e-2 


4.45e-2 


4.89e-2 


5.24e-2 


Med. 


7.61e-5 


5.83e-5 


1.82e-4 


1.69e-2 


1.50e-l 


3.96e-2 


4.41e-2 


6.44e-2 


1.00e-l 


8.12e-2 


U. Qt. 


1.20e-4 


9.35e-5 


2.08e-4 


2.13e-2 


1.95e-l 


4.97e-2 


5.89e-2 


1.17e-l 


1.96e-l 


1.34e-l 


C5 


Mean 


9.41e-5 


5.99e-5 


1.80e-4 


1.65e-2 


1.36e-l 


3.97e-2 


4.63e-2 


7.91e-2 


9.15e-2 


9.96e-2 


St. Dev. 


4.24e-5 


3.13e-5 


6.86e-5 


8.11e-3 


6.83e-2 


1.25e-2 


1.84e-2 


5.48e-2 


6.27e-2 


5.84e-2 


L.Qt. 


5.68e-5 


3.78e-5 


1.31e-4 


1.28e-2 


1.05e-l 


3.05e-2 


3.18e-2 


3.86e-2 


4.91e-2 


5.11e-2 


Med. 


8.96e-5 


5.85e-5 


1.78e-4 


1.69e-2 


1.37e-l 


3.96e-2 


4.59e-2 


5.74e-2 


7.81e-2 


7.24e-2 


U. Qt. 


1.15e-4 


9.01e-5 


2.07e-4 


2.08e-2 


1.67e-l 


5.17e-2 


5.98e-2 


1.04e-l 


9.95e-2 


1.20e-l 



TABLE VII: Summary of Monte Carlo results for parameter estimation errors. We show the mean, standard deviation, median 
and quartiles of the distribution of the error in each parameter. Results are given for a m = 1.4Mq CO inspiralling into a 
M — WOMq IMBH with spin parameter q = 0.3, and assuming four alternative configurations for the detector network, C1-C4, 
as described above. Configuration C5 is the network of three ETs which has been used for all results elsewhere in this paper. 
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Statistics of distribution for error AX in parameter X = 



Model 


ln(m) 


ln(M) 


q 


to 


</>o 


e s 


<t>s 


6k 


4>K 


ln(D) 




Mean 


2.14e-3 


5.01e-4 


2.34e-3 


2.08e-2 


2.34e-l 


1.5332 


2.3089 


2.3442 


2.7840 


2.5105 




St. Dev. 


1.07e-3 


3.51e-4 


1.03e-3 


7.33e-3 


1.30e-l 


2.0017 


3.1318 


1.9584 


2.2266 


1.8864 


CI 


L.Qt. 


1.78e-3 


2.63e-4 


1.95e-3 


1.02e-2 


1.62e-l 


0.1341 


0.2344 


0.8615 


1.1489 


1.1982 




Med. 


2.51e-3 


5.01e-4 


2.34e-3 


1.62e-2 


2.45e-l 


0.6232 


0.8482 


1.9277 


2.0915 


1.8508 




U. Qt. 


3.16e-3 


8.51e-4 


2.63e-3 


2.69e-2 


3.38e-l 


2.1273 


3.1569 


3.4303 


3.8956 


3.2680 




Mean 


1.70e-3 


4.17e-4 


2.19e-3 


1.31e-2 


1.51e-l 


6.08e-2 


9.23e-2 


1.99e-l 


2.28e-l 


2.39e-l 




St. Dev. 


1.15e-3 


3.13e-4 


7.79e-4 


6.46e-3 


7.06e-2 


3.02e-2 


5.75e-2 


1.73e-l 


1.90e-l 


1.86e-l 


C2 


L.Qt. 


1.26e-3 


2.09e-4 


1.82e-3 


8.71e-3 


1.12e-l 


3.50e-2 


4.27e-2 


8.03e-2 


8.81e-2 


9.77e-2 




Med. 


2.29e-3 


3.72e-4 


2.19e-3 


1.28e-2 


1.47e-l 


5.95e-2 


7.01e-2 


1.29e-l 


1.63e-l 


1.63e-l 




U. Qt. 


2.88e-3 


7.24e-4 


2.51e-3 


2.04e-2 


2.08e-l 


7.81e-2 


1.37e-l 


2.51e-l 


2.87e-l 


3.16e-l 




Mean 


1.48e-3 


3.80e-4 


2.09e-3 


1.23e-2 


1.34e-l 


4.38e-2 


5.56e-2 


l.lle-1 


1.31e-l 


1.32e-l 




St. Dev. 


9.07e-4 


2.43e-4 


7.01e-4 


5.76e-3 


4.70e-2 


1.56e-2 


2.52e-2 


9.34e-2 


1.08e-l 


1.02e-l 


C3 


L.Qt. 


1.05e-3 


2.14e-4 


1.78e-3 


8.13e-3 


1.04e-l 


2.99e-2 


3.41e-2 


5.37e-2 


5.60e-2 


6.74e-2 




Med. 


2.14e-3 


3.55e-4 


2.14e-3 


1.25e-2 


1.28e-l 


4.18e-2 


4.95e-2 


7.82e-2 


9.88e-2 


9.70e-2 




U. Qt. 


2.69e-3 


6.76e-4 


2.45e-3 


1.90e-2 


1.73e-l 


5.71e-2 


7.43e-2 


1.35e-l 


1.70e-l 


1.76e-l 




Mean 


1.51e-3 


3.72e-4 


2.09e-3 


1.20e-2 


1.28e-l 


3.97e-2 


5.31e-2 


1.05e-l 


1.21e-l 


1.20e-l 




St. Dev. 


8.06e-4 


2.51e-4 


6.08e-4 


5.09e-3 


4.75e-2 


1.13e-2 


1.96e-2 


8.40e-2 


9.73e-2 


9.51e-2 


C4 


L.Qt. 


1.15e-3 


1.95e-4 


1.70e-3 


7.59e-3 


1.04e-l 


3.16e-2 


3.49e-2 


5.24e-2 


5.71e-2 


6.66e-2 




Med. 


2.09e-3 


3.39e-4 


2.14e-3 


1.23e-2 


1.31e-l 


4.06e-2 


4.85e-2 


8.21e-2 


9.40e-2 


9.82e-2 




U. Qt. 


2.63e-3 


6.56e-4 


2.45e-3 


1.86e-2 


1.62e-l 


4.99e-2 


7.01e-2 


1.39e-l 


1.59e-l 


1.72e-l 




Mean 


1.41e-3 


3.63e-4 


2.09e-3 


1.19e-2 


1.05e-l 


3.95e-2 


5.12e-2 


9.57e-2 


1.08e-l 


1.12e-l 




St. Dev. 


6.39e-4 


1.87e-4 


5.72e-4 


5.07e-3 


5.20e-2 


1.19e-2 


2.03e-2 


6.74e-2 


8.75e-2 


7.67e-2 


C5 


L.Qt. 


1.07e-3 


2.00e-4 


1.65e-3 


7.41e-3 


7.83e-2 


2.96e-2 


3.43e-2 


5.74e-2 


5.75e-2 


5.93e-2 




Med. 


1.58e-3 


3.47e-4 


2.13e-3 


1.24e-2 


9.38e-2 


3.86e-2 


4.79e-2 


6.91e-2 


8.79e-2 


8.82e-2 




U. Qt. 


1.74e-3 


6.03e-4 


2.40e-3 


1.41e-2 


1.22e-l 


4.92e-2 


6.51e-2 


1.14e-l 


1.24e-l 


1.56e-l 



TABLE VIII: As 



Table IVHI but for binary systems with a CO of mass m = 10M© , and a central IMBH of mass M = 5OOM0 . 



Statistics of distribution for error AX in parameter X = 



Model 


ln(m) 


ln(M) 


<7 


to 


00 


0s 


4>s 


9 K 


4>K 


ln(D) 




Mean 


3.90e-4 


6.97e-5 


2.14e-4 


2.75e-2 


3.37e-l 


1.0248 


1.4132 


1.73984 


1.98796 


1.96526 




St. Dev. 


9.60e-5 


1.30e-5 


5.03e-5 


1.12e-2 


2.11e-l 


1.0433 


1.5308 


0.98158 


1.25075 


1.21136 


CI 


L.Qt. 


3.56e-4 


6.24e-5 


1.87e-4 


2.23e-2 


1.33e-l 


0.1554 


0.2162 


0.94480 


0.93815 


0.97202 




Med. 


4.04e-4 


7.02e-5 


2.12e-4 


2.69e-2 


2.76e-l 


0.6771 


0.8102 


1.60984 


1.78872 


1.76174 




U. Qt. 


4.43e-4 


7.79e-5 


2.32e-4 


3.80e-2 


6.02e-l 


1.5840 


2.2538 


2.42813 


2.82516 


2.74045 




Mean 


3.37e-4 


6.92e-5 


2.10e-4 


1.81e-2 


1.38e-l 


5.69e-2 


8.02e-2 


1.93e-l 


2.37e-l 


1.88e-l 




St. Dev. 


1.31e-4 


1.29e-5 


5.00e-5 


9.60e-3 


1.07e-l 


2.85e-2 


5.13e-2 


1.61e-l 


2.05e-l 


1.31e-l 


C2 


L.Qt. 


2.41e-4 


6.23e-5 


1.91e-4 


1.62e-2 


6.26e-2 


3.24e-2 


3.67e-2 


7.81e-2 


9.54e-2 


8.84e-2 




Med. 


3.89e-4 


7.02e-5 


2.11e-4 


1.99e-2 


8.27c-2 


5.02e-2 


6.51e-2 


1.39e-l 


1.54e-l 


1.43e-l 




U. Qt. 


4.23e-4 


7.46e-5 


2.25e-4 


2.57e-2 


1.816-1 


7.23e-2 


1.15e-l 


2.56e-l 


3.16e-l 


2.56e-l 




Mean 


3.23e-4 


6.90e-5 


2.05e-4 


1.62e-2 


1.25e-l 


4.36e-2 


5.96e-2 


1.00e-l 


1.38e-l 


9.97e-2 




St. Dev. 


1.30e-4 


1.18e-5 


4.57e-5 


6.12e-3 


1.06e-l 


1.77c- 2 


3.07e-2 


8.18e-2 


l.lle-1 


6.14e-2 


C3 


L.Qt. 


2.10e-4 


6.23e-5 


1.95e-4 


1.47e-2 


6.30e-2 


2.75e-2 


3.31e-2 


4.84e-2 


6.45e-2 


5.62e-2 




Med. 


3.75e-4 


6.97e-5 


2.06e-4 


1.77e-2 


8.12e-2 


3.96e-2 


5.75e-2 


6.82e-2 


9.11e-2 


7.89e-2 




U. Qt. 


4.13e-4 


7.58e-5 


2.19e-4 


2.23e-2 


2.45e-l 


5.82e-2 


7.43e-2 


1.29e-l 


1.73e-l 


1.39e-l 




Mean 


3.21e-4 


6.85e-5 


2.03e-4 


1.20e-2 


1.23e-l 


4.06e-2 


5.88e-2 


9.12e-2 


1.24e-l 


9.74e-2 




St. Dev. 


1.19e-4 


1.01e-5 


4.43e-5 


5.33e-3 


7.95e-2 


1.47e-2 


2.81e-2 


6.56e-2 


1.03e-l 


6.36e-2 


C4 


L.Qt. 


2.30e-4 


6.25e-5 


1.93e-4 


1.07e-2 


5.88c-2 


2.99e-2 


3.46e-2 


4.64e-2 


6.38e-2 


5.09e-2 




Med. 


3.73e-4 


6.86e-5 


2.04e-4 


1.23e-2 


8.12e-2 


3.96e-2 


5.62e-2 


6.81e-2 


8.71e-2 


7.63e-2 




U. Qt. 


4.13e-4 


7.45e-5 


2.15e-4 


1.73e-2 


2.516-1 


5.17e-2 


6.61e-2 


1.19e-l 


1.66e-l 


1.29e-l 




Mean 


3.17e-4 


6.82e-5 


2.03e-4 


1.14e-2 


9.02e-2 


4.02e-2 


5.75e-2 


8.59e-2 


1.10e-l 


9.68e-2 




St. Dev. 


1.03e-4 


8.51e-6 


4.31e-5 


5.22e-3 


5.32e-2 


1.28e-2 


2.84e-2 


6.22e-2 


9.12e-2 


6.10e-2 


C5 


L.Qt. 


2.21e-4 


6.28e-5 


1.91e-4 


1.02e-2 


7.02e-2 


3.02e-2 


3.23e-2 


4.48e-2 


5.25e-2 


4.91e-2 




Med. 


3.51e-4 


6.80e-5 


2.01e-4 


1.17e-2 


9.34e-2 


4.16e-2 


5.17e-2 


6.14e-2 


8.22e-2 


7.39e-2 




U. Qt. 


4.06e-4 


7.31e-5 


2.11e-4 


1.65e-2 


1.24e-l 


5.09e-2 


6.97e-2 


1.06e-l 


1.50e-l 


1.21e-l 



TABLE IX: As Table IVTT1 but for binary systems with a central IMBH of mass M - 500M Q . 



